Gate-controlled superconductivity in diffusive multiwalled carbon nanotube 
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We have investigated electrical transport in a diffusive multiwalled carbon nanotube contacted 
using superconducting leads made of Al/Ti sandwich structure. We find proximity-induced super- 
conductivity with measured critical currents up to I cm = 1.3 nA, tunable by gate voltage down 
to 10 pA. The supercurrent branch displays a finite zero bias resistance which varies as Ro oc L^m 
with a = 0.74. Using IV-characteristics of junctions with phase diffusion, a good agreement is 
obtained with Josephson coupling energy in the long, diffusive junction model of A.D Zaikin and 
G.F. Zharkov (Sov. J. Low Temp. Phys. 7, 184 (1981)). 



Superconductivity in carbon nanotubes is an intrigu- 
ing subject. To understand it one has to consider many 
facets of modern physics, including Luttingcr liquid be- 
havior owing to strong electron-electron interactions in 
one dimension as well as Kondo physics due to odd, un- 
paired electronic spin . Intrinsic superconductivity has 
been observed in ropes of single walled carbon nanotubes 
(SWNT) @, while proximity induced superconductiv- 
ity has been investigated recently in individual SWNTs 
[H El El The magnitude of the observed supercurrent 
has varied substantially. In Refs. 0, Q respectively, a 
supercurrent on the order of lOx larger and lOx smaller 
than theoretically expected was observed. Morpurgo et 
al. and J0rgenssen et al, on the other hand, did observe 
only increased conductance near zero bias. 

In multiwalled carbon nanotubes (MWNT), supercur- 
rents have been even harder to achieve, presumably due 
to problems with disorder and impurities. Enhanced 
conductance was observed by Buitelaar et al. Q near 
zero bias, which was interpreted in terms of multiple 
Andrcev reflections (MAR) in the presence of inelastic 
processes Q. Recently, proximity induced supercurrent 
has been observed by Kasumov et al. Q as well as by 
Haruyama and coworkers 0, Hoj , most notably in multi- 
shell-contacted tubes grown within nanoporous alumina 
templates [lpj . Here we report proximity-induced super- 
conductivity that is achieved in an individual, diffusive 
MWNT using bulk(side)-contacted samples with Ti/Al 
contacts. Wc find that the supercurrent can be smoothly 
controlled by gate- voltage, via tuning of the diffusion con- 
stant, and a good agreement is obtained using analysis 
based on long, diffusive SNS junctions supplemented with 
phase diffusion effects, modeled in terms of the resistively 
and capacitively shunted junction model (RCSJ). 

Our tube material, provided by the group of S. 
Iijima, was grown using plasma enhanced growth with- 
out any metal catalyst [llj. The tubes were dispersed 
in dichloroethane and, after 15 min of sonication, they 
were deposited onto thermally oxidized, strongly doped 
Si wafers. A tube of 4 /xm in length and 16.6 nm in diam- 
eter was located with respect to alignment markers using 
a FE-SEM Zeiss Supra 40. Subsequently, Ti contacts of 



width 550 nm were made using standard overlay lithog- 
raphy: 10 nm titanium layer in contact with the tube 
was covered by 70 nm Al in order to facilitate proximity 
induced superconductivity in Ti at subkclvin tempera- 
tures. Last, 5 nm of Ti was deposited to prevent the 
Al layer from oxidation. The length of the tube section 
between the contacts was 400 nm. The electrically con- 
ducting body of the silicon substrate was employed as a 
back gate, separated from the sample by 150 nm of SiC>2. 
An AFM image of our sample is displayed in the inset of 
Fig. [U 

On our "dipstick" dilution refrigerator (Nanoway 
PDR50), the samples were mounted inside a tight copper 
enclosure. The measurement leads were filtered using an 
RC filter with time constant of 1 fis at 4.2 K, followed by 
twisted pairs with tight, grounded electrical shields for 
filtering between still and the mixing chamber, while the 
final section was provided by a 0.7-m-long Thcrmocoax 
cable on the sample holder. In the measurements, cur- 
rent bias was employed via a 100 Mfl room-temperature 
resistor when searching for supcrcurrents, while voltage 
bias was employed when looking at multiple Andreev re- 
flections. In the latter case, also the differential resis- 
tance R = dV/dl was recorded using lock-in techniques. 
The gap of the contact material, Ai ead = 139 fJ,eV, was 
found to be suppressed from our regular gap for Alu- 
minum leads (A reg = 200 /icV) by 30 %. Normal state 
results were measured at B = 0.2 T. 

Our nanotube sample is slightly n-type doped ini- 
tially, as in our previous measurements on semiconduct- 
ing MWNTs [12 ■ The back gate capacitance C g = 4 aF 
was deduced from the measured gate period of SET os- 
cillations AV = 40 mV. The total capacitance of the 
4 /xm-long nanotube is estimated to be Cs = 0.4 fF 
which corresponds to a Coulomb energy of Ec ~ 2.5 
K. Since the Fermi-level of the tube shifts according to 
Cg/C^eAVg, the number of channels can be varied only 
in the range of N ~ to a few over the employed gate 
range of Vg = — 10... + 10 V. Most of the observed 
changes in resistance R(V g ) are then due to a variation 
in the distribution of individual transmission coefficients, 
not from the changes in N. There is a part of the change, 
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FIG. 1: (color online) Current as a function of bias voltage V 
at a few values of gate voltage V g = 3.214 V (•); V g = 3.220 

V (o); V g = 3.226 V (■); V g = 3.232 V (□); V g = 3.244 

V (A); V g = 3.344 V (A). The solid straight line displays 
the normal state IV-curve measured in a magnetic field of 
B = 0.2 T at V g = 3.202 V. The arrows A and B illustrate 
the determination of the maximum supercurrent I cm in the 
non-hysteretic and hysteretic cases. The inset on the lower 
right illustrates a magnification of the low voltage part of the 
IV-curve at V g = 3.214 V where clear hysteresis is visible and 
the critical current I cm = 0.15 nA. The inset on the upper 
left displays an AFM image of our sample (scale bar: 1 /im). 
The data were measured in a two-lead configuration on the 
middle section at T — 65 mK. 



up to 50 % at most (dependent on V g ), that can be at- 
tributed to the Kondo effect. Other typical characteris- 
tics for our samples are: the mean free path l m f p ~ 15 
nm and the diffusion constant D ~ vpl m j p = 1 • 10 -2 
m 2 /s, as deduced from the resistance. This means that a 
400 nm long section (see Fig. 1) is expected to be nearly 
in the long SNS junction limit. Wc find several gate volt- 
age values at which the normal state resistance is in the 
range R n = 15 — 20 kf2. The regions with good normal 
state conductance were found to vary by ±20 mV from 
day to day. This is attributed to variations in the back- 
ground charge. Occasionally, we also saw jumps of the 
background charge which were on the order of 0.1 — 0.3 
electrons. 

The measured IV-curves are illustrated in Fig. ^ The 
gate voltage has been varied over 100 mV which corre- 
sponds to a change of 2.5 electrons on the island. The 
shape is seen to change from a state with linear IV- 
characteristics at small bias to a fully "blockaded" one. 
The IV-curves in the first category display a pronounced 
kink, followed by a plateau at voltages in the range of 
1 — 10 /LtV; at intermediate gate voltages even hystere- 
sis is observed. This behavior is identified as proximity 
induced superconductivity with the kink/plateau region 
indicating the maximum measured supercurrent I cm . A 
more detailed view of the small bias regime is given in the 
lower right inset of Fig. ^ There is hysteresis at the inter- 
mediate values of the critical currents. After a maximum 
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FIG. 2: (color online) Maximum supercurrent 1 cm vs. zero 
bias resistance Ro and normal state resistance R n measured 
at T — 80 mK. The filled red triangles and the open green 
circles displays Ro at gate voltage V g = 5.98 . . . 6.02 V and 
3.21 . . . 3.33 V, respectively. The filled blue circles denote R n 
at V 3 = 3.21 . . . 3.33 V. The solid curve is a power law fit with 
I cm cx Rq 1 ' 3 "' while the dashed line represents I cm oc i?^ 4 ' 93 . 
The dotted curve displays the phase-diffusion relation I cm °c 
Ro 1 valid at Ej « k B T. 



hysteresis of Al = I cm — I re tra P = 54 pA at V g — 3.214 
V, the hysteresis is seen to diminish as the current is 
lowered. The maximum product for I cm Rn = 22 /zeV 

"C Ai eal i- 

According to the RCSJ model, there is hysteresis if the 
scaled temperature T = ksT/Ej is small enough and the 
McCumber parameter (3 — (oj p RCtot) 2 >> 1 where ui p is 
the plasma frequency, R is the shunt resistance, and Ctot 
is the total capacitance involved in the plasma oscillation. 
To make an estimate, we take R = Rj = ~ 2 kfl from 
the IV of the junction above the plateau region, Ctot = 

400 fF, and lu p = J~^jf- t , we get an estimate (3 ~ 24 at 

J c o = 5 nA where 7 c0 = '^§ J - is taken at Ej/k B = 120 
mK. This estimate is close to critical damping and no 
hysteresis at T ~ 1 is to be expected Nevertheless, 
a suppression of I c q = 5 nA down to I crn = 1.3 nA takes 
place by thermal fluctuations. When V g is tuned, Rj 
increases more strongly than T, leading to the appearance 
of hysteresis in the IV-curves due to a larger value for 
(3. Eventually the increase of T with lowering Ej takes 
over and the hysteresis disappears with decreasing I cm as 
observed in Fig. Thus, RCSJ model can qualitatively 
explain the main characteristics observed in Fig. ^ 

The supercurrent branch is found to display a finite 
resistance which depends only weakly on bias (see the 
inset of Fig. ^ This weakness of bias dependence distin- 
guishes our results from those of Buitelaar et al and 
of J0rgenssen et al @ who both observed a clear peak 
in the conductance. The dependence of the measured 
supercurrent I cm on Rq = ^-\v=o as well as on the nor- 
mal state resistance R n is given in Fig. [21 I cm could be 
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tuned over two orders of magnitude from 1.3 nA down to 
10 pA, while the normal state resistance increased only 
by a factor of two: from 15 to 35 kfL We find that the 
data can be fitted by a power law behavior I cm oc R^ 135 . 

The relatively large zero bias resistance, Ro > 1-4 kfi, 
is in a gree ment with the ordinary picture of phase dif- 
fusion |l4l fl5| , which may coexist with hysteretic IV- 
characteristics provided that the environment of the junc- 
tion is frequency dependent: at uj p , Z env ~ Zq = 377 
f2 stabilizin g p hase diffusion, while at low frequencies 
Zenv S> Zq |l3j. Ingold et al. have derived for the 
zero bias resistance due to phase diffusion 
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where Iq(x) represents a modified Bcsscl function. Sub- 
sequently, Grabert et al. have shown that Eq. Q is 
rather accurate, within a factor of ~ 2, even when the 
quantum fluctuations are included In our analysis, 
however, we stick to classical phase diffusion because the 
Coulomb energy E c in our samples remains small owing 
to the large, environmental shunting capacitance, which 
yields fcsT / ' E c ~ 50 and negligible corrections from the 
charging effects. 

In Ref. Hfil the IV-characteristics was derived in the 
limit Ej, eV << k B T , according to which there is a sim- 
ple relation between I c q and I cm : 
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Thus, independent of Z env , one expects J, 
overdamped limit. From the maximum value of I cm = 1.3 
nA, we get Ej/ks = 90 mK at T — 65 mK, which is 
at the limit of applicability of Eq. J5J). Eq. (Q yields 
Ro oc EJ 2 in the limit Ej << k B T. Therefore, using 
Eq. @, we get 7 cm oc Rq~ , which is seen to coincide 
quite well with our data in Fig. [21 where this dependence 
is shown as a dotted line. 

The Josephson energy for a long diffusive junction, 
without interaction effects, can be calculated from the 
equation 0,0] 
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3 + 2V2eR r , 



L 



exp 



which is valid in the limit A/exh — > oo when T ~ 
ieTh/ks and where Lt = \JhD jl^ksT . We are em- 
ploying this formula in our analysis as we are not aware 
of any appropriate theoretical treatment of long, inter- 
acting SNS junctions. In SINIS structures, perturbation 
analysis of the interacting case has shown that there are 
loga rithmic corrections that reduce Josephson coupling 
|2(Tj . but this theory is not suitable for our case as the 
contacts have a high transparency. We also neglect res- 
onance effects in our analysis, i.e. the contribution that 
might be related to Kondo effect. 
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FIG. 3: Temperature dependence of I cm (•) and Ro (o) mea- 
sured at V g = 3.489 V. The solid curve is a fit obtained from 
Eqs. © and J3J using D = 2.2 • 1(T 3 m 2 /s, L = 0.4 fim, 
and R„ — 17 kSl. The dashed line displays -Ro calculated 
from Eq. (J using Z en v = 400 Q. and Ej(T) from the above 
T-dependence fit. 




FIG. 4: (color online) Conductance G vs. bias voltage V. 
Gate voltage V g has been stepped over V g = 3.323 . . . 3.355 
V in steps of 2 mV. The dashed vertical lines indicate loca- 
tions for the first Andreev reflection process if governed by 
unrenormalized Ai ea d/e = 139 fieV. 



By combining Eqs. (J2J and ©, we obtain an analyti- 
cal formula for I cm that has only one fitting parameter, 
namely D (or Lt at certain T). In Fig. [21 we compare 
this formula with the measured temperature dependence 
of I cm . The solid curve in Fig. [3 is the result of the 
fit using L T = y/mxnK/T 188 nm. This thermal length 
corresponds to D — 2.2 • 10~ 3 m 2 /s, yielding CTh = 9.0 
fieV |2l|. The discrepancy at the lowest temperatures 
may be an indication that Eq. Q becomes invalid and 
numerical analysis based on the Ivanchenko-Zilberman 
theory should be done. Deviations at the lowest temper- 
atures are also observed between the fit of Eq. JJJ and 
the measured i?o in Fig. 

Fig. 01 displays the differential conductance G for 
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the low conductance IV-curves of Fig. ^ i n the range 
V g = 3.323 . . .3.355 V. A sequence of maxima is seen, 
which are related to multiple Andreev reflections (MAR) 
|23L E3 | . The MAR peaks are more prominent than what 
we would expect for a long, diffusive contact on the ba- 
sis of a recent numerical analysis by Cuevas et al. [22j |. 
The dashed vertical lines indicate locations for the first 
Andreev reflection process if it is governed by unrenor- 
malized Ai eac i/e = 139 fiV. Since Coulomb effects tend 
to shift MAR peaks upwards in voltage, the gap has to 
be modified at the interface by 20 % downwards, roughly 
similar to findings by J0rgensscn et al. Q in SWNTs 
with Al/Ti contacts. Using A = 0.8 ■ Ai ead = 111 /icV 
we find that the main Andreev peaks are located at |A 
and I A, though the latter one does not coincide exactly 
to the expected location at V > 0. 

The weak gate dependence of the MAR lines is quite 
similar to that found in Ref. 0- Between A and 2A there 
is an additional peak that may be connected to Thou- 
less energy. As we approach the supercurrent region by 
decreasing the gate voltage, erh (D) increases and, con- 
sequently, the peak should move towards 2 A j^- In our 
experiment, however, the peak moves towards A. Notice 
also that in the data of Fig. the supercurrent peak near 
zero bias starts to develop before any signs of higher or- 
der Andreev peaks. This seems to contradict with the 
scenario of Vecino et al. Q who argue that inelastic pro- 
cesses enhance conductance due to higher order MAR 
processes and lead to superconductor-like IV-curves with 
hysteresis at small bias. 

Since the number of transmission channels is rather 
small in our sample, it is possible that the subgap trans- 
port is basically dominated by one single channel. This 
might result, especially, from the Kondo resonance that 
is known playa role in the conductance of good qual- 
ity MWNTs pfl. In fact, when comparing the shape of 
the IV curves in Fig. QJwith the calculated IV-curves for 
single-channel S-N-S contacts |2f| rather good agreement 
is obtained for the range of transmissions r = 0.3 — 0.7. 
Therefore, even though a description using the model of 
a long, diffusive junction seems to work well, presumably 
an analysis based on a set of transmission channels would 
yield an even better agreement. 

In summary, we have observed proximity-induced su- 
perconductivity in a diffusive PECVD-grown MWNT. 
The Josephson coupling energy could be tuned by gate 
voltage via a change in the diffusion constant D oc l/R n - 
The model for long diffusive junctions was successfully 
employed for calculating the dependence of Ej on R n 
and T. The measured IV curves (the zero bias resistance 
i?o and maximum supercurrent I cm ) could be understood 
using analysis based on classical phase diffusion, which 
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leads to a decrease of Ir- 



as E 2 j. 



At T = 65 mK, the 



largest obtained Josephson coupling energy was 0.09 K. 
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